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Numerical  Studies  of  the  Physics  and 
Operation  of  LTG  Materials  and  Devices 


Abstract 

This  document  summarizes  Scientific  Research  Associates,  Inc.,  (SRA)  low  tem¬ 
perature  material  studies,  carried  out  under  U.  S.  Air  Force  Office  of  Scientific  Research 
(AFOSR)  Contract  F49620-94-C-0024.  The  study  summarizes  a  model  that  was  devel¬ 
oped  that  is  consistent  with  present  low  temperature  growth  (LTG)  experimental  studies. 
SRA’s  study  included  one-dimensional  transient  simulations  and  two-dimensional  time 
independent  constrained  geometric  studies.  The  broad  aspects  of  the  study  indicate  that 
annealed  LTG  GaAs  is  best  represented  as  material  containing  precipitates  with  charac¬ 
teristics  of  embedded  Schottky  barriers.  These  embedded  barriers  are,  in  turn  surrounded 
by  defects.  Carrier  transport  in  annealed  LTG  GaAs  is  between  the  precipitates,  with  the 
details  determined  by  the  precipitate  spacing,  the  concentration  of  traps,  and  properties  of 
the  surrounding  traps.  The  two  dimensional  studies  provide  numerical  evidence  that  car¬ 
riers  travel  between  precipitates  and  are  influenced,  to  first  order,  by  the  properties  of  the 
surrounding  traps. 
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Numerical  Studies  of  the  Physics  and 
Operation  of  LTG  Materials  and  Devices 

1.  Introduction 

This  document  summarizes  Scientific  Research  Associates,  Inc.,  (SRA)  low  tem¬ 
perature  material  studies,  carried  out  under  U.  S.  Air  Force  Office  of  Scientific  Research 
(AFOSR)  Contract  F49620-94-C-0024. 

Workers  at  SRA  have  been  studying  the  electrical  properties  of  annealed  low 
temperature  grown  GaAs  and  have  developed  a  model  that  is  consistent  with  present  ex¬ 
perimental  studies.  The  broad  aspects  of  the  study  indicate  that  annealed  LTG  GaAs  is 
best  represented  as  material  containing  precipitates  with  characteristics  of  embedded 
Schottky  barriers.  These  embedded  barriers  are,  in  turn  surrounded  by  defects.  Carrier 
transport  in  annealed  LTG  GaAs  is  between  the  precipitates,  with  the  details  determined 
by  the  precipitate  spacing,  the  concentration  of  traps,  and  properties  of  the  surrounding 
traps.  SRA’s  study  included  one-dimensional  transient  simulations  and  two-dimensional 
time  independent  constrained  geometric  studies.  Work  covering  the  first  two  and  one- 
half  years  of  the  study  has  been  published  in  the  Journal  of  Applied  Physics.  This  study 
is  incorporated  as  Appendix  A  to  this  final  report. 

Recent  studies  have  focused  on  the  details  of  the  multidimensional  transport. 
These  studies,  which  are  also  summarized  in  this  document,  provide  numerical  evidence 
that  carriers  travel  between  precipitates  and  are  influenced,  to  first  order,  by  the  proper¬ 
ties  of  the  surrounding  traps.  These  studies  have  gone  beyond  the  embryonic  stage,  but 
additional  work  is  required  to  provide  meaningful  comparison  to  experiment. 

2.  Structure  of  the  Report 

The  report  is  divided  into  five  (5)  sections  including  the  appendix.  Section  3  is  a 
summary  of  the  LTG  studies  including  aspects  of  the  preliminary  two-dimensional  stud¬ 
ies.  Section  4  discusses  recommendations  for  additional  new  multi-dimensional  LTG 
studies.  Appendix  A  is  a  copy  of  the  Journal  of  Applied  Physics  paper,  summarizing  the 
first  two  and  one  half  years  of  the  study. 
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3.  Summary  of  LTG  Studies 

Most  of  SRA’s  LTG  efforts  involved  a  careful  interpretation  of  experiment  in 
light  of  the  simulation  studies  of  annealed  LTG  GaAs.  Both  transients  and  steady  state 
studies  were  performed  with  emphasis  on  one-dimensional  simulations.  A  wide  variety 
of  models  were  examined  to  explain  the  semi-insulating  properties  of  GaAs  LTG  materi¬ 
als.  These  models  included;  (1)  homogeneous  defects;  (2)  clustered  defects,  (3)  Schottky 
precipitates,  (4)  Schottky  precipitates  plus  residual  defects,  (5)  one-dimensional  geome¬ 
tries,  (6)  cylindrical  geometries,  (7)  spherical  geometries,  (8)  transient  photo-excitations, 
and  (9)  transport  through  multiple  precipitates. 

The  broad  conclusion  of  the  study  was  that  the  precipitates  in  GaAs  were  three 
dimensionally  confined  metallic-like  Schottky  barriers,  and  that  these  precipitates  were 
surrounded  by  defects.  The  Schottky  barrier  precipitates  account  for  the  dc  semi- 
insulating  properties  of  LTG  GaAs;  the  transient  properties  require  the  incorporation  of 
traps  in  the  surrounding  material.  The  reasoning  that  led  to  this  conclusion  is  discussed 
in  the  Journal  of  Applied  Physics  paper  contained  in  the  appendix. 

The  major  shortcoming  of  one-dimensional  Schottky  harrier  precipitate  calcula¬ 
tions  is  that  transport  across  the  precipitate  requires  contributions  from  thermionic  emis¬ 
sion  over  the  barrier  or  tunneling  through  the  barrier.  In  actual  devices  carrier  transport 
will  likely  involve  percolation  paths  around  and  between  precipitates.  While  a  suitable 
treatment  of  percolating  paths  requires  three-dimensional  device  simulations,  much  of  the 
essential  physics  can  be  gleaned  from  scaled  two-dimensional  studies.  Such  two- 
dimensional  studies  were  initiated  under  the  current  program  and  are  summarized  below. 
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3.1  Two  Dimensional  Calculations  with  Precipitates  as  Clusters  of  Defects 

For  the  two-dimensional  studies  the  precipitates  were  represented  as  cylinders. 

The  height  of  each  cylinder  extended  to  ±qo  ,  and  so  we  are  interested  in  transport-per- 

unit-device-depth.  To  set  the  stage  for  the  study  we  first  performed  a  number  of  base 

line,  zero  bias  calculations.  For  example,  figure  1  displays  contours  of  constant  density 


◄ -  lOOnm  - ► 


Figure  1.  Density  contours  for  a 
rectangular  array  of  defect 
precipitates  and  a  doping  of 


◄ -  lOOnm  - ► 


Figure  2.  Constant  conduction 
energy  band  contours  for  the  array 
of  figure  1. 


for  a  rectangular  array  of  nine  precipitates,  where  the  precipitates  are  represented  as 
clusters  of  traps.  The  size  of  the  structure  is  lOOnm  on  edge,  with  20nm  separating  the 
centers  of  the  precipitates.  The  precipitates  are  surrounded  by  a  trap  free,  uniform  doping 
concentration  of  lO^Vm^.  As  indicated  in  figure  1,  the  concentration  between  precipitates 
peaks  near  6xlO^Vm^. 

Figure  2  displays  contours  of  constant  conduction  band  energy  corresponding  to 
the  charge  distribution  of  figure  1 .  While  there  is  apparent  structure  in  the  energy  band 
profile,  because  of  the  negligible  depletion  surrounding  the  defect  precipitates,  the  poten¬ 
tial  energy  difference  between  the  center  of  the  precipitate  and  the  surrounding  substrate 
is  a  small  O.leV.  The  situation  for  defect  precipitates  embedded  within  lower  back- 


5 


ground  concentration  regions  results  in  greater  depletion  as  shown  in  figures  3  and  4, 
Here  we  display  contours  of  constant  conduction  band  energy  for  the  same  defect  struc¬ 
ture  as  that  of  figure  1  and  2— the  difference  is  that  the  shallow  donor  density  is  reduced 
to  in  figure  3  and  to  in  figure  4. 


3.2  Two  Dimensional  Calculations  with  Precipitates  as  Schottky  Barriers 

Figures  5  and  6  are  the  Schottky  barrier  equivalents  of  figures  1  and  2.  Here,  we 

have  a  rectangular  array  of  Schottky  barrier  precipitates  each  separated  by  20nm.  The 
barrier  heights  of  the  Schottky  precipitates  are  each  O.SeV.  The  important  physics  differ- 


◄ - ► 

200nm 


Figure  3.  Contours  of  constant  en¬ 
ergy  for  a  rectangular  array  of  de¬ 
fect  precipitates,  surrounded  by  a 
doping  of  10“/m^. 


M -  lOOnm  - ► 


Figure  4.  Contours  of  constant 
enei^  for  a  rectangular  array  of 
defect  precipitates,  surrounded 
by  a  doping  of  10“/m^. 


ence  between  the  results  of  the  defect  precipitates  and  the  Schottky  precipitates  is  that  in 
the  former  the  system  is  conductive,  whereas  for  the  Schottky  precipitate  the  system  is 
still  non-conducting.  (The  origin  of  this  difference  is  discussed  in  detail  in  the  Journal 
of  Applied  Physics  Paper,  see  appendix.)  As  seen  in  figure  5,  the  region  approximately 
halfway  between  the  precipitates  sustains  a  density  is  approximately  eight  orders  of 
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magnitude  less  than  the  shallow  donor  concentration.  In  the  region  surrounding  the  bar¬ 
rier  the  density  is  fourteen  orders  of  magnitude  smaller. 


3.3  Current  Flow  in  2D  Structures:  Precipitates  as  Clusters  of  Defects  and  Why 
Space  Charge  Limited  Current  Appears  in  LTG  Materials 


Figure  6.  Contours  of  constant  con- 

Figure  5.  Contours  of  constant  density  duction  energy  band  for  the  array 

for  a  rectangular  array  of  Schottky  pre-  figure  5. 

cipitates,  surrounded  by  a  doping  of 
lO^W. 


We  were  also  able  to  perform  calculations  that  illustrate  the  current  paths  in  the 
presence  of  precipitates.  These  calculations  are  preliminary  in  that  they  were  performed 
using  the  defect  precipitate  model,  as  the  defect  model  is  a  computationally  more  forgiv¬ 
ing  model.  (Similar  calculations  with  Schottky  barrier  precipitates  await  future  develop¬ 
ment.)  These  calculations,  nevertheless,  shed  light  on  why  the  experimental  characteris¬ 
tics  show  space  charge  limited  properties.  First,  note  that  in  the  Journal  of  Applied 
Physics  article,  see  appendix,  we  speculated  that  transport  in  the  presence  of  precipitates 
would  likely  be  between  the  precipitates  where  the  carriers  would  be  subject  to  the  stan¬ 
dard  bulk  transport  properties-modified  by  the  presence  of  overlapping  depletion  layers. 


For  these  calculations,  whose  results  are  shown  in  figures  7-9,  the  potential  en¬ 
ergy  is  set  to  zero  at  the  top  of  the  structure  and  is  equal  to  +0.5  eV  at  the  bottom  of  the 
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structure.  Figure  7  displays  density  contours  under  conditions  of  applied  bias.  (Figure  8 
displays  the  density  contours  upon  which  we  have  superimposed  the  current  streamlines.) 
First  note  the  distortion  of  the  depletion  layers  under  bias.  There  is  enhanced  depletion 
surrounding  the  layer  closest  to  the  top  of  the  structure.  It  is  in  this  region  where  a  larger 
fraction  of  the  potential  falls,  as  seen  in  figure  9.  Thus  there  would  appear  to  be  a  block¬ 
ing  region  on  the  top  part  of  the  precipitate  structure. 

The  juxtaposition  of  current  streamlines  and  density  contours  in  figure  8  shows 
current  flowing  from  bottom  to  top.  While  we  had  difficulty  in  suitably  registering  the 
current  streamlines  with  respect  to  the  density  contours  the  essential  features  of  transport 
are  clear  from  this  picture.  Current  flow  is  around  the  precipitates. 

Thus  the  broad  picture  of  transport  in  annealed  material  dominated  by  metallic 
precipitates  is  a  follows: 

(1)  Current  paths  are  between  precipitates. 

(2)  Carriers  traveling  between  precipitates  will  sample  the  material  characteristics 
of  the  host  material,  and  be  responsive  to  defects,  traps,  doping  variations,  etc. 
Carriers  traveling  between  precipitates  are  able  to  sustain  bulk-like  trap  filled 
limited  transport.  Transients  will  be  determined  by  the  transport  of  carriers 
between  the  precipitates.  (Recent  work  at  UCSB  (Ibbetson’s  thesis)  suggests 
a  richer  trap  description 

(3)  The  semi-insulating  properties  of  the  material  are  dominated  by  overlapping 
depletion  regions  associated  with  the  Schottky  barrier-like  metallic  precipi¬ 
tates. 
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Density  Contours  Under  Finite  Bias 


Current  Streamlines  and  Density  Contours 


Figure  7.  Constant  density  con¬ 
tours  for  the  parameters  of  figure 
2,  under  conditions  of  current  flow 
at  a  bias  of  0.5  volts. 


Figure  8.  The  constant  density  con¬ 
tours  of  figure  7,  upon  which  current 
streamlines  are  superimposed. 


Figure  9.  Contours  of  constant  en¬ 
ergy  for  the  parameters  of  figure  7. 
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3.4  Dependence  of  2D  Results  on  Precipitate  Configuration 

We  anticipate  that  the  results  of  the  above  discussion  do  not  depend,  at  least  in  a 

qualitative  manner,  on  the  configuration  of  the  precipitates.  To  demonstrate  this  we  be¬ 
gan  a  study  of  the  charge  distribution  associated  with  a  hexagonal  close  packed  array  of 
precipitates.  Figures  10  and  11  illustrate  for  precipitates  represented  as  a  cluster  of  de- 


Figure  10.  Contours  of  constant 
density  for  a  hexagonal  close  packed 
array  of  defect  precipitates  for  a 
doping  of 


Figure  11.  Contours  of  constant  con¬ 
duction  band  energy  for  the  structure  of 
figure  12. 


fects.  The  precipitate  is  similar  to  that  associated  with  figure!,  the  difference  is  the  con¬ 
figuration.  As  can  be  seen  the  distribution  of  potential  and  density  is  qualitatively  similar 
to  that  of  the  rectangular  array.  We  have  also  taken  the  hexagonal  close  packed  configu¬ 
ration  and  subjected  the  structure  to  an  applied  bias.  As  in  the  case  of  the  rectangular  ar¬ 
ray  the  current  path  is  between  the  precipitates.  Thus  qualitatively  similar  behavior  is 
seen  under  bias. 


4.  Recommendations 

Our  picture  of  transport  in  LTG  material  indicates  that  at  a  minimum  a  two  di¬ 
mensional  picture  is  required.  Here  we  have  current  flow  between  the  precipitates,  the 
magnitude  of  which  is  determined  by  the  applied  bias  and  the  overlapping  depletion  re- 
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gions  of  nearest  neighbor  precipitates.  As  seen  by  the  self-consistent  calculations,  repre¬ 
sented  by  figures  7-9,  the  overlapping  depletion  zones  are  modulated  by  the  applied  bias. 
The  picture  is  similar  to  a  rectangular  array  of  PBTs  with  floating  gates,  and  indicates 
that  the  detailed  properties  of  the  material  between  the  precipitates  will  have  a  first  order 
influence  on  transport,  and  subsequent  time  dependent  behavior.  This  picture  requires 
that  we  extend  the  two-dimensional  modeling  to  include  transient  transport  between 
metallic  precipitates. 
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Numerical  studies  of  annealed  non-stoichiometric  low  temperature 
grown  GaAs 

J.  P.  Krestovsky  and  H.  L.  Grubin®’ 

Scienlific  Research  Associates.  Inc.  Glastonbury,  Connecticut  06033 
(Received  23  August  1996;  accepted  19  February  1997) 

The  means  by  which  deep  traps  in  annealed  low  temperature  molecular  beam  epitaxy  grown  GaAs 
are  responsible  for  its  remarkable  semi-insulating  and  short  lifetime  properties  has  been  the  subject 
of  much  discussion.  For  example,  while  low  dc  bias  and  high  speed  phototransient  measurements 
can  be  explained  as  consequences  of  a  homogeneous  trap  distribution  in  the  non-stoichiometric 
material,  experiments  demonstrating  the  presence  of  electrically  active  precipitates  must  be  dealt 
with.  The  study  below,  which  is  numerical,  concludes  that  a  consistent  argument  based  on  the 
presence  of  electrically  active  precipitates,  coupled  to  a  surrounding  distribution  of  traps,  will 
account  for  much  of  the  observed  experimental  phenomena.  O  1997  American  Institute  of Physics. 
[80021-8979(97)01711-8] 


I.  INTRODUCTION 

Over  the  past  several  years  considerable  research  has 
been  performed  in  an  effort  to  understand  the  favorable  in¬ 
sulating  and  optical  properties  of  GaAs  grown  by  molecular 
beam  epitaxy  (MBE)  at  temperatures  in  the  vicinity  of  200 
°C.  Depending  on  the  precise  growth  temperature  and  sub¬ 
sequent  annealing,  low  temperature  (LT)  grown  GaAs  buffer 
layers  have  been  found  to  reduce  backgating,  sidegating  and 
light  sensitivity  in  GaAs  metal-semiconductor  field-effect 
transistors  (MESFETs).**^  Metal-semiconductor-metal 
(MSM)  structures  have  also  been  fabricated  using  LT  GaAs 
and  tested  as  high  speed  photo  conductive  switches.^’"^ 

The  specific  properties  of  LT  GaAs  have  been  associated 
with  the  presence  of  a  very  hi^  density  of  defects,  on  the 
order  of  10^®/cm^  in  the  as-grown  material,  and  the  subse¬ 
quent  precipitation  of  excess  As  after  the  material  is  an¬ 
nealed.  Look  et  cd?  concluded  that  prior  to  the  annealing 
process  the  conductivity  of  the  material  is  due  to  hopping 
conduction  in  a  densely-populated  ‘*EL2”-like  band.  At  an¬ 
neal  temperatures  above  400  ®C,  the  conduction  mechanism 
becomes  one  of  predominantly  frcc-carricr  conduction,  with 
a  corresponding  reduedon  in  the  deep  donor  concentration. 
The  resistivity  of  the  material  rises  from  values  near  14 
fl  cm  for  the  unarmcalcd  material  to  values  in  excess  of 
10^  n  cm  for  armealed  material.  Additional  early  evidence 
of  an  EL2-likc  defect,  related  to  the  Asq.  defect,  was 
reported.^’^  More  recent  electro-optical  measurements  per¬ 
formed  by  Korona  et  ai^  support  this  result  Early  evidence 
of  a  deep  acceptor  with  an  ionization  eneigy  of  £^+03  eV 
was  reported  by  Kaminska  et  al.?  Tliis  is  also  supported  by 
positron  annihilation  evidence  of  deep  acceptors,  as  reported 
by  Haulojirvi  et  al.^^  and  by  the  results  of  Fang  and  Look^^ 
and  McQuaid  et  al.}^ 

The  highly  resistive  nature  of  LT  GaAs  has  been  attrib¬ 
uted  to  the  presence  of  defects  in  the  as-grown  material  and 
the  precipitation  of  arsenic  in  the  annealed  material.  The  fact 
that  the  material  remains  highly  resistive,  even  when  doped 
with  shallow  donors  to  levels  approaching  10*®/cm'’,  lends 


“Electronic  mail:  hal@srai.com 


support  to  the  presence  of  traps  at  densities  approaching 
10*Vcm^  in  the  annealed  material.  The  presence  of  defects 
can  also  explain  the  short  carrier  lifetimes  satisfactorily.  But 
the  material  and  electrical  characteristics  of  uiiannealed  and 
annealed  material  are  different,  and  Warren,  et  ai  postu¬ 
lated  that  As  precipitates  in  the  annealed  material  act,  elec¬ 
trically,  as  buried  Schottky  barriers.  In  this  model  the  deple¬ 
tion  regions  surrounding  adjacent  Schottky  barriers  overlap, 
thus  creating  a  highly  resistive  material.  The  model  of  War¬ 
ren  et  ai  has  received  considerable  support  from  scanning 
tunneling  microscopy  (STM)  studies  of  Feenstra  et  in 
which  LT  GaAs  layers,  grown  at  225  °C,  annealed  at  600 
°C  and  800  °C,  and  n-type  silicon-doped  at  10*’/cm^,  were 
studied.  The  presence  of  a  significant  density  of  midgap 
states  induced  by  the  precipitates  demonstrated  that  the  pre¬ 
cipitates  were  capable  of  pinning  the  Fermi  level,  and  the 
observation  of  local  depiction  sunounding  the  precipitates 
indicated  that  the  precipitates  were  charged. 

It  has  also  been  proposed*^’*^  that  the  precipitates  could 
act  as  recombination  centers,  thus  accounting  for  the  short 
carrier  lifetimes.  Look^^  indicates  that  for  precipitate  densi¬ 
ties  near  1 X  10*Vcm^  the  lifetime  would  be  about  0.8  ps. 
This  suggestion  has  been  countered  by  the  observation  that 
the  carrier  lifetimes  in  the  unarmcalcd  material,  in  which 
there  are  no  precipitates,  are  of  the  same  order  in  the  an¬ 
nealed  material.  Whitakcr^^  has  reported  lifetimes  on  the  or¬ 
der  of  180-260  fe  for  samples  grown  at  190  and  200  ®C, 
with  and  without  aimealing  at  600  ®C  (for  10-15  min).  Be¬ 
cause  the  lifetimes  are  inversely  proportional  to  the  defect 
densities,  Whitaker’s  results^’  also  suggest  that  annealing  re¬ 
sults  primarily  in  the  precipitation  of  the  As^  defects  which 
form  the  precipitates,  while  As^  and  Vq*  defect  densities 
should  remain  relatively  unaltered..  However,  other  data  sug¬ 
gests  that  As^  and  Vo*  densities  decline  when  the  material 
is  annealed. 

One  feature  is  clear  from  the  above  discussion:  traps  are 
responsible  for  the  behavior  of  non-stochiometric  GaAs.  For 
annealed  material  there  is  the  additional  feature  of  precipi¬ 
tates.  What  are  the  properties  of  these  precipitates?  How  may 
they  be  modeled?  How  do  they  affect  transport?  The  pres- 


7326 


J.  Appl.  Phys.  81  (11).  1  June  1997 


002 1 -8979/97/81  (1 1  )/7326/1 9/$1 0.00 


©  1997  American  Institute  of  Physics 


cncc  of  precipitates  suggests  that  current  paths  in  LT  GaAs 
will  have  percolation  properties,  and  that  simple  bulk  resis¬ 
tivity  formulas  may  not  provide  an  adequate  explanation  of 
device  measurement.  Indeed,  the  fact  that  several  apparently 
different  models  for  the  material  properties  of  LT  GaAs  can 
explain  such  features  as  the  low  field  resistivity  measure¬ 
ments  and  the  carrier  lifetime  measurements,  serves  to  high¬ 
light  the  limitation  of  explaining  a  rich  set  of  material  prop¬ 
erties  on  the  basis  of  current-voltage  and  photo-optical 
transient  measurements.  The  results  of  Feenstra  et  al.  also 
force  us  to  recognize  the  possibility  that  geometrical  effects 
are  important. 

Numerical  simulation  is  an  ideal  tool  to  examine  these 
problems,  with  a  full  scale  study  requiring  realistic  device 
lengths,  transients  and  three  dimensions.  Such  studies  await 
future  efforts.  For  the  present  stuefy,  implementing  numerical 
simulation  procedures  to  study  both  LT  GaAs  material  and 
device  characteristics,  we  have  analyzed  specific,  but  limited 
device  and  material  phenomena.  The  simulation  results  fol¬ 
low  a  pattern  leading  to  the  conclusion  that  spherical  (three- 
dimensional  confined)  metallic-like  precipitates,  coupled  to  a 
surrounding  distribution  of  traps,  will  account  for  much  of 
the  observed  experimental  phenomena. 

II,  THE  MODEL 

The  model  is  based  upon  the  semiconductor  drift  and 
diffusion  equations,  including  trap  kinetics.  The  simulation  is 
based  upon  the  numerical  solution  of  these  equations.  We 
have  used  the  model  to  explore  the  electrical  properties  of 
devices  with  and  without  precipitates.  The  effects  of  precipi¬ 
tates  were  studied  initially  in  a  one-dimensional  planar  ge¬ 
ometry,  and  then  extended  to  single  precipitates  in  spherical 
geometry.  Specifically,  we  have: 

(1)  examined  the  kinetics  of  traps  in  materials  with  ho¬ 
mogeneous  and  inhomogeneous  trap  distributions, 

(2)  modeled  precipitates  as  a  confined  high  density  of 
traps;  this  model  is  referred  to  as  the  defect  precipitate  model 
(DPM), 

(3)  modeled  precipitates  as  Schottlty  barriers;  this  model 
is  referred  to  as  the  Schottky  barrier  precipitate  model 
(SPM), 

(4)  examined  an  enhanced  DPM  which  incorporates  de¬ 
fects  in  the  surrounding  material,  and 

(5)  examined  an  enhanced  SPM,  vriiich  also  incorporates 
defects  in  the  surrounding  material. 

Throu^  the  use  of  numerical  simulation  we  were  able  to 
isolate  the  consequences  of  each  model.  A  “roadmap”  of  the 
stucty,  displayed  as  Table  I,  serves  as  a  guide  to  the  many 
different  types  of  simulations  performed.  The  “yes”  and 
“no”  comments  identify  the  ability  of  the  model  to  explain 
experimental  observations. 

III.  THE  EQUATIONS  AND  THE  MODEL  FORMULATION 

The  simulations  of  LT  GaAsiAs  are  based  on  the  drift 
and  diffusion  equations,  Poisson’s  equation,  and  rate  equa¬ 
tions  for  deep  level  donors  and  acceptors.  The  continuity 
equations,  using  the  drift  and  diffusion  relations  for  current 
densities,  are: 


TABLE  L  Sumnuuy  of  results  of  the  study.  The  yes*  md  *iio*  comments 
identify  the  ability  of  the  model  to  explain  experimental  observations. 


Model  studied 

Precipitates 

Resistivity 

Phototransients 

Uniform  trap 
distribution 

No 

Yes 

Yes 

Defect  precipitate 

Yes 

Yes,  marginal 

No 

model 

geometrical 

constraints 

Schottky  precipitate 
model 

Yes 

Yes 

No 

Enhanced  DPM 

Yes 

Yes,  marginal 
geometrical 

Yes 

constraints 

Enhanced  SPM 

Yes 

Yes 

Yes 
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(0 


dP 
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with  the  constitutive  equations: 

i„=-e{ft„NV<t>„-D„VN),  (3) 

Jp  =  +  e(/ip/’V^^+DpVP).  (4) 

In  equations  (1)  and  (2),  N  and  P  are  the  electron  and  hole 
densities.  In  equations  (3)  and  (4)  and  Pp  are,  respec¬ 
tively,  the  electron  and  hole  mobilities  (with  values  given  in 
Appendix  B)  and  D„  and  Dp ,  the  respective  diffusivities.  G 
and  R  represent  generation  and  recombination,  which  will  be 
presented  subsequently.  For  homogeneous  material  with 
variations  in  the  doping  concentration  the  potentials  in  equa¬ 
tions  (3)  and  (4)  differ  by  the  value  of  the  band  gap.  For  the 
embedded  Schottky  barrier  model  the  precipitates  arc  repre¬ 
sented  as  material  different  fiiom  the  host  material  (vath  a 
band  structure  consistent  with  the  Anderson  model)  and  the 
potentials,  4>„  and  <f>p ,  are  given  by: 

(5) 

€  € 


X  kgT 

4>p^il'+J+- - —rnN„.  (6) 

Here  is  the  electron  affinity.  Eg  the  bandgap,  and  Me  and 
are  the  density  of  states  for  die  conduction  and  valence 
bands.  The  quantity  ^  is  obtained  fiem  Poisson’s  equation: 

V .  (ev  =e[(N-No -N;)-(P+Po+NJ)].  (7) 

where  Nq  and  Pq  are  the  densities  of  shallow  dopants,  and 
and  are  ionized  deep  donors  and  acceptors.  The  total 
number  of  deep  donor  and  deep  acceptor  traps  are  given  by: 

(8) 

(9) 

In  the  following  discussion  we  deal  with  both  linear  and 
spherical  geometries.  For  spherical  geometries,  the  diver¬ 
gence  in  equation  (1)  is:  ( l/r^)(^/<?r)(/'^),  and  the  left  hand 
side  of  equation  (7)  becomes:  {\/r^){d/  dr){r^  €){dl  dr). 
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The  kinetics  of  capture  and  emission  fix)m  donors  and 
acceptors  determine  the  recombination  rates  for  electrons 
and  holes  introduced  into  equations  (1)  and  (2). 

G„-R„  =  C„{N^,N:-NN^^)  +  C:[NiN%-NN*] 

+  G„  (externaf) ,  (10) 

=  - PN^HCpiPiN*  - PN%) 

+  G p{exiernal)  ^  (1 1) 

where  C„  and  denote  capture  rate  of  electrons  from  ac- 
ceptors  and  donors,  respectively,  and  C’  and  Cp  denote  cap¬ 
ture  rates  of  holes  from  acceptors  and  donors,  respectively. 
The  superscripts  +,  and  on  the  donor  and  acceptor  den¬ 
sities  indicate  an  ionized  or  neutral  state.  Wc  have  also  con¬ 
sidered  band-to-band  recombination  through  B{pn—n])^ 

with  ^=10“^°cmVs.  Band-to-band  recombination  effects 

are  negligible,  except  for  the  very  long  time  scales,  and  are 
ignored  in  most  of  the  calculations. 

The  terms  A^i  and  P  \  are  obtained  from  equilibrium  con¬ 
ditions,  and  under  non-degenerate  conditions  are  approxi¬ 
mated  by  the  equations: 


N\=N^  cx^-[{E^-E^)lknT], 

Ari=N,  exp-[{E,-E^)/kBTl 

(12) 

Pl=N^  exp-[(£„-£:j/*:fir]. 

(13) 

where  is  the  conduction  band  energy,  the  valence 
band,  the  energy  level  of  deep  donors,  and  E^  that  of  the 
deep  acceptors. 

The  transients  associated  with  the  filling  and  emptying 
of  the  donor  and  acceptor  states  are  governed  by  rate  equa¬ 
tions  different  from  those  of  equations  (1)  and  (2).  For  a 
single  acceptor  level  and  a  single  donor  level  these  equations 
are: 

(14) 

dN^ 

-^=C:{NiNl-NN*)+Cp(PNl-PiN*).  (15) 

The  models  studied  using  the  above  equations  are: 

A.  The  homogeneous  model: 

The  material  sustains  a  uniform  distribution  of  deep 
level  donors  0.69  cV  below  the  conduction  band  and  deep 
level  acceptors  03  eV  above  the  valence  band.  Concentra¬ 
tion  dependence  was  studied,  as  well  as  the  effect  of  shallow 
dopants.  Uniformly  distributed  donors  and  acceptors  are 
most  appropriate  for  unannealed  material,  but  are  included  in 
this  discussion  of  annealed  material. 

B.  The  defect-precipitate  model  (DPM): 

This  model  assumes  that  the  precipitate  sites  in  annealed 
LT  GaAs  differ  from  the  surrounding  material  only  through 
an  extraordinarily  high  local  concentration  of  defects.  Since 
the  material  constituting  the  precipitates  is  different  from  the 
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host  material,  the  conclusions  of  the  study,  with  this  assump¬ 
tion,  are  regarded  as  only  qualitatively  relevant.  Assuming 
precipitates  on  the  order  of  60  A  in  diameter  at  a  density  of 
lO^Vcm^,  and  a  total  defect  density  of  10‘Vcm^,  the  density 
of  defects  within  a  precipitate  would  be  8.85X  lO*Vcm^  if  all 
the  defects  were  concentrated  within  the  precipitates  (see 
Appendix  B).  For  the  simulations  we  assume  that  these  pre¬ 
cipitates  are  periodically  distributed  throughout  the  material 
at  a  spacing  of  200  A,  which  is  consistent  with  a  density  of 
10*^  precipitates/cm^.  (The  dependence  of  the  results  on  pre¬ 
cipitate  density  is  also  discussed.) 

C.  Computational  Schottky  barrier  precipitate  model: 

Most  device  studies  that  incorporate  Schottky  barriers, 
c-g,,  field-effect  transistor  (FET)  studies,  do  not  deal  whh  the 
specific  properties  of  the  Schottky  barrier;  rather,  they  treat 
the  effect  of  the  barrier  as  a  boundary  condition.  However, 
this  is  not  the  case  when  the  object  of  study  is  the  Schottky 
barrier.  For  the  present  study,  while  we  are  not  interested  in 
the  detailed  properties  of  the  Schottky  barrier,  treating  the 
Warren  et  ai  Schottky  barrier  model  as  a  boundary  condi¬ 
tion  would  not  be  fhiitful.  Instead,  when  simulating  precipi¬ 
tates  as  buried  Schottky  barriers,  the  precipitate  was  repre¬ 
sented  as:  (1)  a  zero  band  gap  semiconductor  with  a  specified 
offset  voltage  with  respect  to  GaAs,  in  which  (2)  the  electron 
and  hole  density  are  equal  and  high  (lO^^/cm^)  and  where 
(3)  everywhere  within  the  material,  except  at  the  region  sur¬ 
rounding  the  boundary  of  the  zero  band  gap  material  and 
the  GaAs  material,  there  was  local  charge  neutrality.  In  ad¬ 
dition,  by  setting  the  hole  mobility  and  diffrisivity  to  zero 
within  the  precipitate,  the  holes  become  immobile  and  be¬ 
have  like  ionized  sites.  The  electron  mobility  and  diffrisivity 
are  set  to  high  values,  which  3dclds  the  desired  low  resis¬ 
tance,  A  generic  representation  of  this  model  is  illustrated  in 
Figure  1,  with  the  Schottky  region  to  the  left  of  the  origin. 
While  there  are  no  traps  within  the  modeled  Schottky  barrier 
region,  the  model  permits  a  thin  layer  of  deep  level  donors 
and  acceptors  at  the  semiconductor-Schottky  interfece  to  rep¬ 
resent  surfecc  recombination  effects.  Because  the  model  is 
new,  a  series  of  representative  calculations  arc  discussed  in 
Appendix  A. 

D.  Enhanced  defect  precipitate  model: 

This  model  starts  with  the  defect  precipitate  model  and 
adds  traps  in  the  region  between  the  precipitates.  These  traps 
may  have  spatially  varying  properties. 

E.  Enhanced  Schottky  precipitate  model: 

This  model  starts  with  the  Schottky  precipitate  model 
and  adds  traps  in  the  region  between  the  precipitates.  These 
traps  may  have  properties  that  vary  spatially. 

IV.  SPATIALLY  CONSTANT  RESULTS:  RESISTIVITY 
AND  THE  FERMI  LEVEL  IN  THE  ABSENCE  OF 
PRECIPITATES 

While  annealed  LTG  GaAs  contains  precipitates  and  the 
measured  resistivity  is  a  property  of  percolating  paths  as  well 
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FIG.  1.  Generic  representation  of  a  zero-band  gap  metal  semiconductor 
configuration.  The  electron  affinity  varies  discontinuously  across  the  inter¬ 
face. 


FIG.  2.  For  the  case  where  the  dependence  on  Fenni  energy 

(solid  line)  and  resistivity  (dashed  line)  on  the  concentration  of  shallow 
donois  (lower  abscissa).  For  the  case  of  zero  shallow  donois,  the  depen¬ 
dence  of  Fermi  energy  (solid  triangles)  as  function  of  NyN\  (upper  ab¬ 
scissa). 


as  the  defect  kinetics,  the  fact  that  simple  models  can  ac¬ 
count  for  some  of  the  low  field  properties  of  non- 
stochiometric  GaAs  is  addressed  in  this  section.  The  resistiv¬ 
ity  of  the  material  was  computed  for  a  uniform  distribution 
of  traps.  Guided  by  various  experimental  investigations,  we 
studied  material  with  donors  0.69  cV  below  the  conduction 
band  and  acceptors  03  cV  above  the  valence  band.  The  re¬ 
sistivity  of  the  homogeneous  material  is  H(Nefji„"^Pe^p) 
where  the  densities  are  obtained  fiom  the  Idnctic  equations 
subject  to  a  condition  of  diarge  neutrality.  The  electron  and 
hole  densities  in  equilibrium  arc  determined  by  the  Fermi 
level,  which  in  turn  in  diis  study  is  dominated  by  the  position 
of  the  donor  trap,  the  ratio  of  donor  and  acceptor  traps,  and, 
for  shallow  donors,  the  ratio  of  shallow  donors  to  acceptor 
traps.  This  is  displayed  in  Figure  2. 

Figure  2  contains  three  plots.  The  solid  line  displays  the 
relative  position  of  the  Fermi  energy  as  a  function  of  the 
ratio  of  shallow  donors  to  total  acceptors  (lower  abscissa)  for 
the  case  vs^ere  the  ratio  of  deep  donors  to  deep  acceptors  is 
N^/N^  =  30.  This  calculation  was  undertaken  because  of  the 
experimental  observation  that  annealed  LT  GaAs  remains 
semi-insulating,  even  in  the  presence  of  a  high  concentration 
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of  shallow  donors.  The  solid  line  results  are  approximately 
accounted  for  by  the  expression: 

Ec-EF=Ec-EB+kgT/n{l-No/N^) -kgT/n{NyNl) 

(16) 

and  show  that  for  shallow  donor  densities  less  than  O.IA/^  the 
conduedon  band  remains  very  close  to  the  undoped  equilib¬ 
rium  position  of 0.6028  cV  above  the  Fermi  level.  To  under¬ 
stand  this  result,  we  note  that  in  the  absence  of  shallow  do¬ 
nors  the  deep  acceptors  are  completely  ionized,  receiving  an 
electron  fiom  the  deep  donors:  Np  —  .  As  shallow  donors 

arc  introduced  (they  are  assumed  to  be  completely  ionized 
and  degeneracy  is  not  considered)  the  ionized  deep  donors 
are  initially  neutralized  by  an  amount  equal  to  the  shallow 
doping  density.  For  shallow  donor  densities  below  O.W^, 
the  concentration  of  ionized  deep  donors  is  sufficient  for  the 
material  to  remain  hi^ly  resistive;  for  shallow  donor  levels 
in  excess  of  the  Fermi  level  moves  up  toward  the 

conduction  band. 

The  relative  insensitivity  of  the  Fermi  energy  to  the  con¬ 
centration  of  shallow  donors  is  reflected  in  values  of  the 
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resistivity,  also  shown  in  Figure  2  by  a  dashed  line.  For  these 
calculations,  where  we  assumed  a  mobility  of  8000 
cm^A^  s  for  electrons  and  400  cm^A^  s  for  holes,  the  resistiv¬ 
ity  exceeds  lO’  Hem.  The  results  shown  in  Figure  2  are 
typical  of  the  results  obtained  for  other  ratios  of  donor  to 
acceptor  trap  levels.  The  primary  difference  is  the  location  of 
Fermi  energy  for  low  shallow  donor  doping. 

For  later  discussion  we  include  a  plot  of  the  Fermi  level 
as  a  function  of  the  ratio  of  donor  traps  to  acceptor  traps  for 
the  case  of  zero  shallow  donors.  This  is  displayed  by  the 
closed  triangles  (upper  abscissa)  in  Figure  2. 

The  above  results  show  that  a  homogeneous  distribution 
of  traps  could  account  for  the  dc  semi-insulating  properties 
of  LT  GaAs:As.  From  the  experimental  data  (ignoring  pre¬ 
cipitates),  after  anneal  the  level  of  deep  donor  defects  is  on 
the  order  of  10**/cm^,  and  fiem  Figure  2  we  can  conclude 
that  the  material  would  retain  the  semi-insulating  properties 
up  to  a  shallow  donor  level  approaching  10*Vcm^.  For  cases 
where  the  ratio  of  deep  donors  to  acceptor  is  closer  to  l.O, 
the  shallow  donor  level  will  be  even  higher  before  the  mate¬ 
rial  becomes  conducting. 

V.  DC  CHARACTERISTICS  FOR  NIN  STRUCTURES 
WITH  A  HOMOGENEOUS  DISTRIBUTION  OF 
TRAPS 

The  first  level  of  complication  is  the  A^-doped- 
insulator-A^-doped  (NIN)  structure  with  band  conduction. 
We  point  out  an  even  greater  complication  proposed  by  Luo 
et  where  conduction  in  annealed  LT  GaAs  is  a  result 

of  hopping  between  precipitates  (we  will  discuss  this  later). 
Historically,  one  of  the  earliest  numerical  studies  of  band 
conducting  semi-insulating  properties  of  NIN  GaAs  was  that 
of  Horio  et  who  concentrated  on  compensated  material 
with  an  region  20  fim  long.  A  more  recent  study  of 
relevance  to  the  present  article  is  that  of  Ibbetson  et  cd.^ 
with  a  2-/an-long  annealed  /  LT-GaAs  region.  The  Horio 
et  cd}^  structure,  which  is  significantly  longer  than  that  for 
which  non-stochiomctric  GaAs  is  presently  grown,  is  ger¬ 
mane  to  the  present  discussion  because  it  highlights  the  sig¬ 
nificance  of  trap-filled  limited  conduction  as  a  possible  trans¬ 
port  mechanism. 

For  the  structure  of  Horio  et  and  fior  an  /  Ic^r 
containing  deep  donors  and  deep  acceptors  at  a  ratio, 
=30,  as  in  Figure  2,  we  obtained  current  versus  volt¬ 
age  (see  Figure  3).  (Calculations  with  parameters  appropriate 
to  Ibbetson  et  al}^  are  displayed  in  Figure  4.)  Two  sets  of 
calculations  were  performed  on  the  structure  of  Horio 
et  identified  by  symbols:  A^J=6X10'Vcm^  and 

6X10^Vcm^.  No  shallow  acceptors  are  in  the  calculation. 
The  contact  regions  arc  doped  with  shallow  donors  to  a  level 
of  2X  10^/cm^.  For  comparison  we  repeated  Horio  et 
(Figure  3,  lines)  for  EL2  deep  donor  traps  and  shallow  ac¬ 
ceptors,  at  a  fixed  ratio  of  1  Lll;  simulations  were 

performed  for  =  5  X  10*^/cm^,  2X10'Vcm^  and 

5X10'Vcm\ 

While  the  details  may  differ,  all  five  curves  in  Figure  3 
display  qualitatively  similar  behavior.  In  the  low  bias  region, 
electrons  injected  into  the  /  region  are  trapped  and  the  resis¬ 
tivity  of  the  structure  is  determined  primarily  by  the  resist!  v- 


Bias.  log,o(V) 


FIG.  3.  (a)  Structure  for  the  NIN  calculations,  (b)  Current-voltage  relation 
with  as  a  parameter.  Calculations  identified  by  Horio  tt  at.  (see  Ref.  20) 
arc  for  a  ratio  of  donor  traps  to  shallow  acceptors  of  11.11.  Calculations 
represented  by  symbols  are  for  a  ratio  of  donor  traps  to  acceptor  traps  of  30. 
For  all  calculations  the  capture  coefficients  of  Horio  et  al.  (see  Ref  20), 
were  used. 

ity  of  the  /  region.  Until  the  trap-filled  limit  is  approached, 
the  resistivity  remains  relatively  constant  As  the  trap-filled 
limit  is  approached,  the  traps  can  accept  only  a  significantly 
reduced  number  of  additional  electrons,  and  a  small  increase 
in  bias  results  in  a  flooding  of  the  I  layer  with  electrons  and 
a  corresponding  dramatic  drop  in  resistivity.  Consequently, 
the  current  increases  rapidly,  with  small  bias  increases.  Be¬ 
yond  the  trap-filled  limit  the  velocity  saturates  and  the  cur¬ 
rent  is  governed  by  the  space  charge  in  the  /  region,  which 
continues  to  increase  with  bias.  The  details  follow. 

For  the  curves  labeled  Horio  et  al.:  at 
=  5  X  lO^Vcm^,  due  to  the  low  trap  density,  veiy  few 
electrons  need  be  injected  into  the  /  region  before  all  the 
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FIG.  4.  Current- voltage  characteristics  for  a  2  /xm  NIN  structure  with  pa¬ 
rameters  similar  to  that  of  Ibbctson  et  al.  (see  Ref.  21).  Results  are  also 
shown  for  an  /  region,  3  /xm  long 


traps  are  occupied;  the  resistivity  is  small  There  is  a  soft 
transition  between  the  ohmic  and  space  charge  limited  re¬ 
gions.  When  the  trap  density  is  increased  to  Np-2 
XlO*Vcm^,  there  is  substantially  hi^er  resistance  in  the 
low  bias  region,  and  the  low  bias  and  space  charge  limited 
regions  are  delineated  by  a  sudden  increase  in  current  at  the 
trap-filled  limit  For  a  trap  density  of  Nq— 5  X  10*Vcm^  the 
resistance  is  further  increased;  because  additional  bias  is 
needed  to  provide  enough  electrons  to  neutralize  the  traps, 
the  presence  of  velocity  saturation  begins  to  affect  the  IV 
curve  in  the  form  of  saturation  below  the  trap-filled  limit 

The  N^N^=30  calculations  with  acceptor  traps  (sym¬ 
bols  in  Figure  3)  are  qualitatively  similar  to  those  of  Horio 
et  oL?^  The  resistivity  increases  as  the  concentration  of  do¬ 
nor  traps  increases.  The  effects  of  velocity  saturation,  prior 
to  trap-filled  limited  transport  at  higher  donor  trap  concen¬ 
trations,  delays  the  reduction  in  resistivity  associated  with 
trap-filled  limited  transport  Differences  with  Horio  et 
e.g,,  comparing  the  latter  resistance  for  N^=5  X  lO^Vcm^  to 
the  deep  acceptor  trap  calculations  with  AA^=6XlO'Vcm^, 
are  attributed  to  differences  in  the  Fermi  level. 

The  study  of  Ibbctson  et  cd?^  was  an  experimental  and 
analytical  investigation  of  space  charge  limited  currents  in 
nonstoichiometric  GaAs.  The  /  region  was  2  fim  thick  and 
was  separated  from  an  buffer  by  a  20  A  AlAs  layer.  The 
experimental  results  of  Ibbctson  et  show  apparent  trzp- 
filled  limited  operation.  The  qualifier  arises  because  there 
was  neither  a  discussion  of  breakdown,  nor  of  the  length 
dependence  of  the  voltage  at  which  the  current  shows  the 


sharp  increase,  Luo  etai^^  argue  that  the  threshold  voltage 
for  trap-filled  transport  varies  as  the  square  of  the  length  of 
the  /  region,  whereas  breakdown  varies  linearly  with  length. 
Ibbetson  et  using  a  set  of  boundary  conditions  that 
were  equivalent  to  pure  drift  current  far  from  the  junction, 
computed  /-  V  curves  that  displayed  saturation,  but  none  of 
the  trap-filled  limit  conductivity  properties  seen  in  Figure  3, 

We  performed  an  NIN  calculation  for  the  Ibbetson 
et  ai^^  structure  (see  Figure  4).  In  our  study  the  AlAs  layer 
was  ignored,  the  /  region  was  2  pxn  long,  and  the  two  end 
N**^  regions  were  each  0.5  pm  long,  with  a  density  of 
10*Vcm^.  For  this  calculation  the  donor  traps  were  situated 
al  0.73  eV  below  the  conduction  band,  the  mobility  of  the 
electrons  was  1900  cm^/V  s,  =  10**/cm^,  A^^  =  2.4 
X  l0*^/cm^,  which  is  an  order  of  magnitude  smaller  than  that 
used  in  the  simulation  study  of  Ibbetson  et  ai?^  The  accep¬ 
tors  were  treated  as  traps  with  densities  chosen  to  provide 
resistance  values  similar  to  experimental  values.  Addition¬ 
ally,  velocity  saturation  in  our  calculation  is  8  X  10^  cm/s  as 
opposed  to  5  X  10^  cm/s  for  Ibbetson  et  ai?^  These  param¬ 
eters  yield  the  results  shown  in  Figure  4.  (The  calculations 
were  repeated  with  shallow  acceptors;  there  was  no  change . 
in  the  results.) 

The  Figure  4  calculations  display  low  bias  resistivity  re¬ 
sults  similar  to  those  of  Ibbetson  et  ai?^  The  calculations  are 
classic  trap-filled  limit  characteristics,  with  the  high  current 
levels  at  a  voltage  threshold  somewhat  higher  than  that  of  the 
experimental  results.  But  the  results  of  Figure  3  suggest  that 
suitable  tweaking  of  the  input  trap  densities  could  yield 
matching  voltage  thresholds.  The  length  dependence  of  the 
threshold  for  the  onset  of  trap-filled  current  was  also  exam¬ 
ined.  The  longer  the  /  region,  the  higher  the  value  of  the 
threshold  voltage  for  high  current  values.  This  is  seen  in  the 
current— voltage  curve  for  a  structure  with  a  3  pm  long  / 
layer.  For  calculations  wth  /  layer  lengths  of  2, 3,  and  4  ftm, 
the  length  dependcce  of  the  voltage  at  which  the  current 
increased  catastrophically  is  given  approximately  by 
which  is  somewhat  less  than  the  classic  space 
charge  limited  b^avior  which  varies  as  the  square  of  the 
device  length.  The  role  of  velocity  saturation  in  the  depen¬ 
dence  of  the  threshold  voltage  must  be  assessed. 

While  the  simple  model  without  precipitates  appears  to 
account  for  the  broad  IV  characteristics  of  NIN  LT  GaAs 
structures,  it  is  premature  to  conclude  that  transport  in  LT 
GaAs  is  trap  filled  and  then  at  higher  bias  levels  space  charge 
limited  transport  Important  unknowns  here  are  the  capture 
coefficients.  The  calculations  of  Figure  4  were  performed 
using  the  capture  coefficients  of  Horio  et  al}^  (see  Appendix 
B  for  values).  We  then  repeated  the  calculations  of  IV  with 
changes  in  the  capture  coefficients  (O.l  X  standard,  lOX  stan¬ 
dard,  lOOX  standard),  with  the  resulting  simulations  showing 
significant  changes  in  the  hole  density.  But  the  magnitudes 
of  the  hole  density  were  small  and  did  not  change  the  input 
into  Poisson’s  equation,  nor  did  they  contribute  to  the  total 
current  in  any  significant  manner.  As  a  result,  the  IV  char- 
acteristics  were  unchanged  Additionally,  for  a  given  set  of 
deep  donor  trap  levels  the  results  were  insensitive  to  the 
properties  of  the  acceptor  traps.  For  the  default  set  of  capture 
coefficients,  the  IV  was  unchanged  when  the  deep  acceptor 
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FIG.  5.  Transient  response  of  electrons  and  holes  in  the  center  of  an  NIN 
structure  following  the  switch-on  of  100  mW/cm^  of  band  gap  radiation,  and 
the  subsequent  swiich-off  transient. 

traps  were  replaced  by  shallow  acceptor  traps.  This  latter 
result  is  likely  to  be  altered  in  the  presence  of  shallow  donors 
of  any  consequent  density.  If  there  is  arching  that  this  pre¬ 
liminary  model  suggests,  it  is  that  dc  measurements  are  not 
sensitive  to  some  of  the  material  properties  of  annealed  non- 
stoichiometric  GoAs,  Transients  should,  of  course,  be  sensi¬ 
tive  to  capture  coefficients,  and  we  deal  with  this  next. 

VL  THE  HOMOGENEOUS  DEFECT  MODEL; 
TRANSIENT  RESPONSE  FOLLOWING 
PHOTOGENERATION 

As  seen  below,  numerical  simulation  teaches  that  the 
transient  response  of  LT  GaAs  to  photoexcitation  is  sensitive 
to  material  properties  We  first  illustrate  this  for  the  NIN 
structure  with  the  dimensions  of  Figure  3,  with  the  /  layer 
modeled  as  LT  GaAs  with:  (i)  a  homogeneous  distribution  of 
defects,  (ii)  capture  coefficient  values  of  Horio  et  al.^^  (iii) 
nI/N^=2  and  (iv)  ^^o  =  10^Vcm^  Zero  bias  was  main¬ 
tained  across  the  structure,  and  the  /  layer  was  subjected  to  a 
uniform  generation  rate  equivalent  to  100  mW/cm^  of  band 
gap  radiation.  The  transient  response  of  the  electron  and  bole 
density  in  the  center  of  the  device  was  monitored.  Figure  5 
shows  the  electron  and  hole  transients  for  the  swheh-on  and 
switch-off  transient.  As  can  be  seen,  after  approximately  1  ps 
the  electron  density  has  reached  a  steady  density  level  at  a 
value  approximately  three  orders  of  magnitude  greater  than 
that  without  illumination.  For  times  less  than  100  ps  holes 
closely  follow  the  electrons,  and  the  densities  of  both  are 
governed  by  the  generation-recombination  times  of  the  con¬ 
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FIG.  6.  As  in  Figure  5  but  for  the  ionized  donors  and  acceptors. 


tinuity  equations  (1)  and  (2).  The  trap  rate  equations  have 
longer  time  scales  than  those  of  electrons  and  holes,  and 
prior  to  100  ps  have  not  responded  to  the  higher  levels  of 
electrons  and  holes.  The  response  of  the  deep  donors  and 
acceptors  is  shown  in  Figure  6.  The  electron  density  relaxes 
back  to  its  initial  value  about  I  ps  efter  the  illumination  is 
switched  off.  The  electron  lifetime  appears  to  be  of  the  order 
of  I  ps. 

Beyond  100  ps,  the  ionized  donor  population  begins  to 
decrease  slightly  in  response  to  the  excess  electrons.  Since 
the  kinetics  of  the  trap  rate  equations  for  the  specified  energy 
levels  and  capture  coefficients  dictate  that  Nq  ,  the  ion¬ 
ized  acceptor  concentration  also  declines.^  As  a  result,  the 
excess  optically  generated  holes  can  no  longer  recombine  at 
the  same  rate  as  the  electron  and  the  hole  density  begins  to 
rise  again.  A  new  equilibrium  is  finally  reached  after  about  1 
ms,  where  the  hole  concentration  is  approximately  2.5 
X  10^/cm^.  This  value  is  the  difference  between  the  ionized 
acceptors  and  donors,  i.e.,  0=NQ—Njf+h—e 
where  e<h<N^  . 

During  switch-offi  the  level  of  ionized  traps  does  not 
begin  to  return  to  the  dark  levels  until  almost  a  microsecond 
has  elapsed.  But  the  most  dramatic  feature  of  this  figure  is 
that  the  transients  for  the  deep  donors  and  acceptors,  as 
shown  in  Figure  6,  are  almost  identical  (to  the  resolution  of 
the  figure  they  appear  exactly  the  same).  This  is  because  of 
the  relatively  large  levels  of  the  ionized  traps,  compared  to 
the  density  of  the  photogenerated  carriers,  particularly  the 
holes.  From  Figure  5  recall  that  during  the  long  time  scale 
portion  of  the  hole  transient  the  hole  density  increases  from 
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FIG.  7.  As  in  Figure  5  with  the  following  change:  the  hole  capture  cocfTi- 
cieni  was  increased  to  the  value  of  the  electron  capture  coefficient. 


FIG.  8.  As  in  Figure  5,  with  the  following  changes;  the  hole  capture  coef¬ 
ficient  was  increased  to  the  value  of  the  electron  capture  coefficient  and  the 
energy  levels  of  the  donor  and  acceptor  traps  were  placed  equally  from  the 
conduction  and  valence  bands. 


lO^/cm^  to  lO^Vcm*’,  which  represents  a  small  perturbation 
compared  to  the  ionized  defect  densities,  >^ich  are  of  the 
order  of  10*^/cm^.  Thus,  the  levels  of  the  ionized  defects 
change  only  slightly  from  the  equilibrium  levels  and,  due  to 
the  long  time  scales  over  which  the  changes  occur,  they  re¬ 
main  in  near  equality  with  each  other.  The  small  overall 
change  in  the  ionized  defect  density  which  docs  occur  is 
made  in  an  effort  to  maintain  charge  neutrality. 

The  origin  of  the  secondary  time  scale  observed  in  Fig¬ 
ure  5  was  traced  to  the  value  of  the  hole  capture  coefficient, 
Cp .  By  increasing  this  coefiScient  to  a  value  equal  to  that  of 
the  electron  coefficient,  C„,  tiic  results  shown  in  Figure  7 
were  obtained,  where  we  still  observe  a  secondary  hole  time 
scale;  however,  it  docs  not  result  in  a  significant  increase  in 
the  hole  density.  The  decay  of  the  holes,  however,  is  still 
significantly  longer  than  friat  of  the  electrons.  Another  con¬ 
tribution  to  the  differences  between  the  electron  and  hole 
lifetime  was  found  to  be  due  to  the  relative  position  of  the 
deep  level  traps  in  the  cnagy  band.  By  placing  the  deep 
level  donors  and  acceptors  equally  from  the  conduction  and 
valence  bands,  respectively,  the  results  of  Figure  8  were  ob¬ 
tained.  Under  these  conditions,  both  electron  and  hole  life¬ 
times  are  close  to  1  ps,  and  no  secondary  time  constants  are 
apparent. 

How  do  the  above  results  translate  into  an  understanding 
of  the  LT  devices  of  the  type  discussed  by  Ibbetson  et  ali^^ 
The  results,  of  course  suggest  that  transient  measurements 
are  necessary  to  fingerprint  non-stoichiometric  materials.  To 
illustrate  this,  an  academic  problem  was  simulated:  the  tran¬ 


sient  response  of  the  2  /zm  long  structure  of  Figure  4  to 
optical  radiation.  The  results  of  this  calculation,  shown  in 
Figure  9,  display  the  sensitivity  of  the  switch-off  transient  to 
the  value  of  the  capture  coefficients.  Individual  calculations 
are  for  (i)  normal  capture  coefficients,  (ii)  coefficients  an 
order  of  magnitude  larger,  and  (iH)  coefficients  two  orders  of 
magnitude  larger.  For  the  case  of  the  nominal  capture  coef¬ 
ficients  the  band  gap  radiation  is  at  a  level  of  100 
mW/cm^.  The  radiation  was  adjusted  for  variations  in  the 
capture  coefficients  to  enable  the  initial  densities  to  be  the 
same.  The  calculation  of  Figure  9  demonstrates  that  tran¬ 
sients  are  sensitive  to  the  capture  coefficients,  even  vben  the 
IV  curves  are  not  The  transients  are  also  sensitive  to  the 
density  of  tiic  traps,  even  when  the  IV  curves  are  not  (sec 
Figure  10).  At  the  very  least,  these  calculations  suggest  that 
dc  measurements  are  not  sufficient  to  characterize  non- 
stoichiometric  GaAs. 

While  the  above  transient  calculations  demonstrate  the 
sensitivity  of  device  performance  to  the  time  constants,  the 
consequences  of  differences  in  the  electron  and  hole  time 
constants  will  be  a  significant  factor  in  assessing  device  per¬ 
formance,  To  illustrate  this  a  simple  friought  experiment  was 
conducted:  For  an  NIN  structure  under  bias  we  examined  the 
response  to  a  photoexcitation  of  10  ps  with  a  repetition  rate 
of  10  ps.  Over  ten  cycles  of  10  mW  photoradiation,  the  elec¬ 
trons  appear  to  relax  to  a  value  close  to  their  initial  value. 
However,  the  hole  density  does  not  relax  to  the  initial  state  at 
the  beginning  of  each  cycle,  due  to  the  longer  relaxation 
time.  As  a  result,  over  a  large  number  of  cycles,  the  excess 


J.  Appl.  Phys.,  Vol.  81,  No.  11.  1  June  1997 


J.  P.  Kreskovsky  and  H.  L.  Grubin  7333 


FIG.  9.  Transient  response  of  electrons  and  holes  in  the  center  of  the  2  /mi 
/  region  of  the  NfN  structure  of  Figure  4  following  the  switch-off  of  band 
gap  radiation.  Capture  cross  section  was  varied  (lOx  denotes  an  order  of 
magnitude  increase  in  cross  section,  etc.).  The  radiation  intensity  was  varied 
such  that  the  initial  carrier  densities  were  the  same. 

holes  accumulate  and  the  current  transient  exhibits  drift. 
When  the  hole  lifetime  is  the  same  order  as  the  electrons, 
both  electrons  and  holes  respond  equivalently  to  the  light 
pulse  and  the  current  transients  are  repetitive.  These  results 
suggest,  at  the  very  least,  that  the  pulse  duration  and  the 
repetition  rate  must  be  adjusted  for  the  slowest  carrier. 

VIL  TWO-DIMENSIONAL  CONTRIBUTIONS 

While  one-dimensional  simulations  are  useful,  two- 
dimensional  effects  are  expected  to  play  an  important  role  in 
the  design  of  devices.  Such  a  two-dimensional  study  would 
cany  us  beyond  the  scope  of  this  discussion,  but  some  com¬ 
ments  are  in  order.  For  example,  when  considering  precipi¬ 
tates,  multidimensional  transport  would  result  in  carrier  flow 
between  precipitates  which,  depending  upon  the  level  of  in¬ 
jection,  would  begin  to  reduce  the  size  of  any  depletion  zone 
surrounding  the  precipitate.  For  the  homogeneous  distribu¬ 
tion  of  traps  considered  here,  the  two-dimensional  effects  are 
concerned  with  the  distribution  of  photogenerated  carriers.  In 
the  one-dimensional  studies,  the  photogeneration  is  assumed 
to  be  uniform  in  a  direction  parallel  to  the  surface.  In  two 
dimensions  flie  generation  decays  with  distance  from  the  sur¬ 
face  due  to  absorption  of  photons.  In  preliminary  studies, 
with  the  photon  flux  and  generation  terms  given  as,  e.g., 
r  =  and  G  =  ar  where  To  is  the  photon  flux  trans¬ 

mitted  through  the  surface  (Fq  is  computed  from  the  optical 
power  =  where  v=c/X,  a  is  the  optical  absorp¬ 


FIG.  10.  Transient  response  of  electrons  and  holes  in  the  center  of  the  2  /im 
/  region  of  the  NIN  structure  of  Figure  4  following  the  switch-off  of  100 
mW/cm^  band  gap  radiation.  Parameter  variation  is  the  acceptor  trap  den¬ 
sity. 


tion  coefficient,  and  x  is  the  distance  of  penetration  measured 
along  the  optical  path,  in  this  case  normal  to  the  device  sur¬ 
face)  the  optically  generated  charge  density  shows  a  similar 
exponential  variation.  However,  at  low  power  levels  the 
charge  density  was  not  great  enough  to  affect  the  potential 
distribution  across  the  device.  Thus,  while  there  is  an  x  di¬ 
rection  in  optically  generated  charge,  the  resulting  current  is 
directed  only  in  the  direction  normal  to  the  contacts.  At 
higher  optical  power  levels,  the  density  of  the  photogener¬ 
ated  charge  would  increase  and  at  some  point  would  affect 
the  potential  distribution,  creating  an  x  direction  field.  The 
current  paths  under  such  conditions  will  be  more  complex 
than  those  addressed  here. 


VIIL  INVESTIGATION  OF  PRECIPITATE  MODELS: 
DEFECT  VERSUS  BURIED  SCHOTTKY 

While  the  previous  discussion  demonstrates  that  ignor¬ 
ing  the  presence  of  precipitates  and  assuming  a  homoge¬ 
neous  distribution  of  donor  and  acceptor  defects  at  high  con¬ 
centration  can  apparently  account  for  such  characteristics  as 
resistivity,  transients,  and  shallow  donor  dependence,  the 
presence  of  precipitates  must  be  included.  We  consider  the 
two  different  precipitate  models  mentioned  earlier:  the 
defect-precipitate  model,  and  the  Schottky  barrier  precipitate 
model. 
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FIG.  11.  Onc-dimcnsional  planar  structure  for  initial  precipitate  simulations. 


IX.  DEFECT-PRECIPITATE  MODEL  (DPM),  PLANAR 
GEOMETRY 

The  defect-precipitate  model  assumes  that  the  precipitate 
sites  in  annealed  LT  GaAs  contain  high  local  concentrations 
of  defects.  Assuming  precipitates  on  the  order  of  60  A  in 
diameter  at  a  density  of  10**^/cm^,  and  a  total  defect  density 
of  10*Vcm\  the  density  of  defects  within  a  precipitate  would 
be  8.85  X  lO'Vcm^  if  all  the  defects  were  concentrated  within 
the  precipitates.  We  further  assume  that  these  precipitates  are 
periodically  distributed  throughout  the  material  at  a  spacing 
of  200  A,  which  is  consistent  with  a  density  of  10*^ 
precipitates/cm^ . 

Initial  simulations  were  performed  using  a  one- 
dimensional  planar  geometry  to  illustrate  key  qualitative  fea¬ 
tures  of  the  problem.  Two  precipitates  were  modeled,  as 
shown  in  Figure  1 1.  The  ratio  of  donor  to  acceptor  defects  in 
the  precipitate  regions  was  set  to  yield  an  equilibrium  Fermi 
level  0.8  eV  below  the  conduction  band,  which  ftom  Figure 
2  requires  approximately  the  same  density  for  each  type  of 
trap.  The  electron  distributions  for  shallow  doping  densities 
between  2  X  lO^Vcm^  and  2  X  lO'^/cm^  are  shown  in  Figure 
12.  These  results  indicate  that  the  inter-precipitate  material 
would  remain  semi-insulating  until  the  doping  level  reached 
2  X  lO'^/cm^,  assuming  that  all  shallow  donors  were  ionized. 
We  have  not  included  the  effects  of  degeneracy  here.  The 
clusters  of  defects  making  up  the  precipitate  are  able  to  trap 
the  majority  of  charge  from  the  shallow  doped  region  pro¬ 
vided  the  shallow  doping  level  does  not  exceed  a  value  that 
depends  upon  the  trap  density  in  the  precipitate.  At  a  doping 
level  of  1 X  lO^^/cm^  the  peak  electron  density  is  on  the  or¬ 
der  of  only  lO^Vcm^. 

Figure  13  shows  the  level  of  ionized  donors  and  accep¬ 
tors  in  the  precipitate  regions  for  a  shallow  donor  level  of 
2X10*’/cm^,  at  which  the  structure  should  be  regarded  as 
conductive.  To  interpret  this  figure,  it  should  be  understood 
that  under  undoped  conditions  the  densities  of  ionized  donor 
and  acceptor  traps  in  the  precipitates  are  the  same,  and  equal 
to,  the  total  acceptor  trap  density  (acceptor  trap  density  is 
below  that  of  the  donor  trap  density).  Thus,  the  figure  shows 
that  the  deep  donors  near  the  surface  of  the  precipitate  are 
being  neutralized  by  trapping  an  electron  from  the  shallow 


FIG.  12.  Electron  distribution  between  and  within  the  vicinity  of  defect 
precipitates.  For  this  calculation  E^— E cV,  A^p  =  8.84x  IO*^/cm^. 

1.014,  and  the  concentration  of  shallow  donors  was  varied. 

donor.  Continued  increase  in  the  shallow  donor  level  results 
in  further  neutralization  of  the  ionized  deep  donors.  The  re- 
suits  displayed  in  Figures  12  and  13  also  depend  upon  the 
inter-precipitate  separation.  Additional  calculations  indicate 
that  as  this  separation  increases  the  peak  doping  density  at 
which  the  material  remains  resistive  is  reduced.  While  not 
shown,  the  potential  energy  distribution  accompanying  the 
above  charge  distribution  is  similar  to  that  at  a  Schottky 
barrier.  For  DPM,  the  effective  barrier  height  is  dependent 
on  the  density  of  the  shallow  dopants  surrounding  the  clus- 
since  the  trap  kinetics  within  the  cluster  determine  the 
electron  density  within  the  cluster,  vriiidi  in  turn  depends 
upon  the  surrounding  doping  level. 

These  above  results  suggest  that  precipitates  consisting 
of  a  high  concentration  of  defects  can  also  account  for  the  dc 
sexni-insuladng  properties  of  LT  GaAs  (see  Table  I).  But 
these  planar  defect  precipitate  models  can  only  be  viewed 
as  qualitative,  since  the  assumption  of  a  planar  geometry  is 
erroneous.  Planar  geometries  allow  excess  depiction  (com¬ 
pared,  e.g.,  to  spherical  geometries)  because  the  number  of 
free  electrons  available  to  be  trapped  and  the  number  of  deep 
level  trapping  sites  vary  linearly  with  distance  from  the 
interface  between  the  precipitate  and  the  inter-precipitate 
material.  Spherical  geometries,  discussed  below,  provide  a 
more  realistic  representation  of  a  single  precipitate,  since  the 
number  of  electrons  to  be  depleted  would  increase  with  the 
cube  of  the  distance  fiom  the  precipitate  boundary,  and  the 
number  of  trap  sites  would  decrease  with  the  cube  of  the 
radius  as  we  move  towards  the  center  of  the  precipitate. 
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FIG.  13.  Distribution  of  ionized  defects  for  the  Figure  12  calculation,  and 
;^o  =  2XlO'’/cni^ 


Thus,  spherical  geometiy  would  not  result  in  such  high  lev¬ 
els  of  depletion  for  the  same  density  of  defects  in  the  pre¬ 
cipitate. 

X.  SCHOTTKY  BARRIER  PRECIPITATE  MODEL  (SPM) 

We  repeated  the  simulations  of  the  planar  DPM  with  the 
precipitate  modeled  as  a  Schott!^  barrier.  The  modeling  of 
the  Schottky  barrier  was  previously  discussed,  with  several 
illustrations  provided  in  Appendix  A.  Figure  14  shows  the 
distribution  of  free  electrons  between  two  Schottky  barriers 
for  the  same  shallow  doping  levels  used  in  the  DPM  simu¬ 
lations.  The  barrier  height  was  taken  as  0.8  eV.  The  results 
are  remarkably  similar  to  those  of  the  DPM.  At  a  shallow 
doping  level  of  1 X  l0*’/cm^  the  peak  electron  densi^  be¬ 
tween  the  barriers  is  also  about  1 X 10*  Vcm^.  However,  at  a 
doping  level  of  2X  10*’/cm^  the  Schottky  barrier  precipitate 
model  (SPM)  yields  a  peak  density  of  just  over  1 
X  10*Vcm^,  whereas  the  DPM  yields  a  significantly  higher 
peak  density,  close  to  I  X  l0*^/cm^.  Charge  neutrality  is 
maintained  in  this  structure  by  the  accumulation  of  electrons 
on  the  “metal”  side  of  the  “metal”-semiconductor  interface 
in  the  SPM,  whereas  in  the  DPM  ionized  donor  defects  were 
neutralized,  exposing  ionized  acceptor  defects.  (Note:  in  a 
true  metal,  whose  density  is  much  higher  than  that  repre¬ 
sented  here,  the  amount  of  charge  accumulation  on  the  metal 
side,  while  consistent  with  charge  neutrality,  would  be  neg¬ 
ligible  compared  to  the  nominal  metallic  density.) 


FIG.  14.  Electron  distribution  between  and  within  the  vicinity  of  Schottky 
barrier  precipitates.  For  this  calculation  the  barrier  height  was  0.8  cV,  and 
the  concentration  of  shallow  donors  was  varied. 


XI.  SPHERICAL  MODEUNG  OF  PRECIPITATES: 
SCHOTTKY  BARRIER  VERSUS  DEFECT  MODELS 

In  view  of  the  qualitative  similarity  found  in  planar  ge¬ 
ometry  simulations  of  the  two  precipitate  models  and  be¬ 
cause  planar  geometries  are  unrealistic,  we  turned  our  efforts 
to  modeling  of  precipitates  using  spherical  coordinates.  As 
was  stated  before,  the  DPM  is  able  to  deplete  the  adjacent 
semiconductor  material  based  on  the  trap  sites.  In  the  planar 
configuration  the  charge  depleted  increases  lineariy  and  the 
trap  sites  available  increase  lineariy  with  distance  fiom  the 
interface  between  the  precipitate  and  the  semiconductor  ma¬ 
terial  In  spherical  geometry,  the  charge  depleted  decreases 
with  the  cube  of  the  depiction  region  width  and  die  trap  sites 
decrease  in  a  similar  manner.  Using  SPM,  the  charge  de¬ 
pleted  is  profxjrtional  to  the  surface  area  of  the  Schottky 
barrier  and  we  expect  to  sec  a  greater  deviation  between  the 
two  models.  However,  while  it  was  possible  to  investigate 
two  adjacent  precipitates  in  a  planar  configuration  and  verify 
that  depletion  regions  would  overlap,  the  spherical  geometry 
used  here  (we  are  not  performing  three-dimensional  simula¬ 
tions)  allows  us  to  mode!  only  a  single  precipitate,  then, 
based  on  planar  studies,  extrapolate  the  results  to  spherically 
confined  precipitates. 

XII.  SPHERICAL  PRECIPITATE-DEFECT  SIMULATION 

To  show  the  geometric  effects  we  first  computed  the 
charge  and  potential  energy  distributions  for  a  single  precipi¬ 
tate  60  A  in  diameter,  with  N^p^lN\  —  2  and  the  donor  trap 
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density  in  the  precipitate  set  at  8,84X  This  yielded 

a  Fermi  level  near  midgap  (£<--£^^=0.69  eV)  under  equi¬ 
librium  conditions.  The  material  surrounding  the  precipitate 
was  doped  with  shallow  donors  at  a  level  of  1  X  lO’Vcm^. 
Comparisons  of  the  spherical  computations  and  planar  ap¬ 
proximations  are  shown  in  Figures  1 5a  and  15b  on  linear  and 
log  scales.  For  the  spherical  geometry  the  DPM  is  unable  to 
deplete  the  surrounding  material.  Within  15  A  of  the  precipi¬ 
tate  surface  (a  radius  of  45  A  from  the  center)  the  density  of 
electrons  in  the  surrounding  material  is  nominally  1 
XlO’Vcm^  At  a  radius  of  100  A  the  electron  density  is 
6.75X  10*Vcm^.  This  suggests  that  if  precipitates  were 
spaced  200  A  apart,  the  minimum  electron  density  resulting 
from  the  overlapping  of  depletion  regions  would  be  on  the 
order  of  3  X  lO^Vcm^  or  greater  for  a  shallow  doping  level  of 
lX10*Vcm^,  This  is  drastically  different  from  the  planar 
case,  where  significant  depletion  is  shown  to  extend  well 
beyond  100  A.  Calculations  were  also  performed  for  a  shal¬ 
low  donor  density  of  1  X  10*’/cm^.  While  the  depletion  re¬ 
gion  extends  further  than  for  the  1  X  lO^Vcm^  case,  at  100  A 
the  density  is  still  on  the  order  of  3  X  10*^/ cm^,  which  would 
render  the  material  conductive.  When  the  trap  ratio  is  in¬ 
creased  to  1.014,  with  the  same  donor  concentra¬ 

tion  of  8.84X  10‘Vcm^,  the  Fermi  energy  increases  to  0,8  cV 
and  the  depletion  region  broadens,  as  seen  in  Figure  15b.^^ 
However,  even  at  this  deeper  Fermi  level,  the  electron  den¬ 
sity  at  100  A  radius  is  over  4.3  X  10*^/cm^  when  the  sur¬ 
rounding  material  is  doped  to  1  X  lO’Vcm^.  This  would  be 
very  close  to  being  on  the  verge  of  the  boundary  between 
conducting  and  semi-insulating  material  when  the  interaction 
of  several  precipitates  is  considered-  From  this  we  can  con¬ 
clude  that  if  precipitates  of  defects  are  responsible  for  the 
semi-insulating  properties  of  LT  GaAszAs,  then  the  semi- 
insulating  properties  would  be  critically  dependent  on  the 
concentration  of  defects  within  the  precipitate. 

XIII.  SPHERICAL  SCHOTTKY  BARRIER  SIMULATIONS 

Simulations  for  the  SPM  in  spherical  coordinates  show 
qualitatively  similar  density  profiles  to  the  DPM,  as  shown  in 
Figure  16.  However,  there  are  important  differences  between 
the  two  models.  In  the  DPM,  the  effective  barrier  height  is 
set  by  the  position  of  the  equilibrium  Fermi  level,  which 
varies  as  die  traps  in  the  precipitate  are  filled.  In  the  SPM, 
the  barrier  height  remains  fixed  regardless  of  the  doping  con¬ 
dition.  The  SPM  is  thus  fer  more  linear  and,  as  a  result,  we 
expect  greater  depletion  than  with  the  DPM.  From  Figure  16 
we  observe  that  at  a  shallow  doping  density  of  l0**/cm^  the 
free  electron  concentration  at  a  radius  of  100  A  is  about 
3.4X  10*^/cm^,  whereas  for  the  defect  model  with  the  equi¬ 
librium  Fermi  level  at  0.8  eV,  the  electron  density  was  4.3 
X  10^^/cm^  at  the  same  doping  level.  This  is  over  an  order  of 
magnitude  greater.  It  is  thus  apparent  that  the  SPM  would 
yield  a  more  resistive  material  than  the  DPM,  unless  signifi¬ 
cantly  higher  trap  levels  were  present  in  the  DPM  precipi¬ 
tate. 

The  DPM  simulations  show  that  the  depletion  region  is 
very  sensitive  to  the  Fermi  level  in  the  precipitate.  SPM 
simulations  show  an  almost  negligible  change  when  the  bar- 


Distance  from  center  of 
precipitate  (Angstroms) 


FIG.  15.  Comparison  of  spherical  and  planar  precipitate  depletion  re¬ 
gions  using  the  DPM  with  a  shallow  donor  concentration  of  =  10**/ cm^ 
and  a  deep  trap  concentration  of  A/^  =  8.84X  lO'^/cm^.  (a)  Linear  scale, 
(b)  logarithmic  scale.  For  the  solid  and  dotted  lines  £'(--£^/r=0.69  eV, 
for  the  dashed  lines  £^—^^=0-8  eV. 
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FIG.  16.  Comparison  of  electron  density  distribution  for  spherical  and  pla¬ 
nar  depletion  regions  using  SPM.  For  the  solid  and  dotted  lines  the  barrier 
height  is  0.8  eV.  For  the  dashed  line  the  barrier  height  is  0.7  eV. 


rier  height  is  changed,  (see  Figure  16,  dashed  line).  In  Figure 
16,  the  density  of  electrons  at  100  A  is  3,8  X 10*^/ cm^  for  an 
0.8  cV  barrier  and  5.4X  10*^/ cm^  for  the  0.7  eV  barrier.  The 
difference  between  the  DPM  and  SPM  results  is  as  follows: 
when  undoped,  the  position  of  the  Fermi  level  is  approxi¬ 
mately  the  same  in  both  models.  However,  in  the  DPM  when 
the  material  is  doped,  the  migration  of  electrons  into  the 
precipitate  and  neutralization  of  ionized  donor  sites  results  in 
a  movement  of  the  Fermi  level  towards  the  conduction  band, 
reducing  the  depletion  width.  This  feature  is  absent  in  the 
SPM. 

The  level  of  traps  required  in  the  DPM  to  yield  results 
similar  to  the  SPM  is  shown  in  Figure  17,  where  it  is  appar¬ 
ent  that  even  with  trap  densities  3  times  greater  than  would 
seem  reasonable,  the  DPM  still  falls  short  of  the  Schottky 
barrier  result  at  barrier  heights  of  0.7  and  0.8  eV. 

XIV,  PHOTO-OPTICAL  RESPONSE  OF  SPM 

With  reference  to  Tabic  I,  three  models  have  been  ex¬ 
amined.  It  was  demonstrated  that: 

(1)  a  model  in  which  a  homogeneous  distribution  of 
traps  is  introduced  can  account  for  both  the  resistive  and 
optical  properties  of  LT  GaAs; 

(2)  a  planar  geometry  model  of  bodi  the  DPM  and  SPM 
can  also  account  for  the  semi-insulating  properties  of  LT 
GaAs; 


Distance  from  center  of 
precipitate  (Angstroms) 


FIG.  17.  Electron  density  distribution  for  DPM  and  SPM.  with  difTcrent 
barrier  heights  and  different  trap  densities. 


(3)  in  a  model  with  spherical  geometry  precipitates  the 
DPM  framework  is  very  sensitive  to  the  density  of  traps,  and 
is  at  the  limit  of  explaining  the  semi-insulating  properties  of 
LT  GaAs; 

(4)  a  spherical  SPM  model  can  easily  account  for  the  dc 
semi-insulating  characteristics  of  the  material. 

Based  upon  these  features,  a  study  of  the  optical  re¬ 
sponse  of  LT  GaAs  within  the  SPM  framework  was  under¬ 
taken  (similar  calculations  with  the  DPM  framework  produce 
the  same  results).  This  investigation  was  undertaken  using  a 
single  precipitate.  Optical  generation  was  introduced  and  the 
switch-on  transient  state  was  obtained;  the  optical  generation 
was  then  turned  off  and  the  transient  monitored.  The  electron 
density  100  A  from  the  precipitate  center  was  observed  as  a 
function  of  time.  For  the  case  where  the  materia!  surround¬ 
ing  the  precipitate  was  trap  free,  it  was  immediately  appar¬ 
ent  that  the  optically  generated  carriers  could  only  be  col¬ 
lected  by  ambipolar  diffusion  through  the  boundary  of  the 
computational  domain,  and  that  the  time  scales  would  be 
considerably  longer  than  those  seen  experimentally.  Band- 
to-band  time  scales  also  fall  into  the  long  time  scale  cat¬ 
egory.  In  view  of  the  long  time  constants,  we  felt  that  the 
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FIG.  18.  Transient  response  of  electron  density  at  a  point  100  A  from  the 
center  of  a  Schottky  barrier  precipitate,  as  a  function  of  residual  trap  density 
surrounding  the  precipitate.  Optical  power  was  adjusted  to  yield  constant 
initial  electron  density. 

simple  SPM  was  likely  to  be  inadequate,  and  that  enhance¬ 
ments  to  the  SPM  and  DPM  were  necessary,  in  the  form  of 
residual  traps  in  the  material  surrounding  the  precipitate, 

XV.  PHOTO-OPTICAL  RESPONSE  AND  CARRIER 
LIFETIME  FOR  ENHANCED  DPM  AND  SPM 

A  series  of  simulations  were  performed  for  the  SPM 
with  residual  traps  in  the  surrounding  material.  In  the  sur¬ 
rounding  material  AAj/iV^=2  and  the  residual  donor  trap 
density  varied  from  (no  residual  traps)  to  1  X  10*Vcm^.  The 
surrounding  materied  was  void  of  shallow  donors,  (For  com¬ 
panion  DPM  simulations,  the  defect  level  in  the  precipitates 
was  also  A^^/AaJ=2,  with  the  donor  traps  set  at  8.84 
X  lO^Vcm^.)  The  optical  illumination  was  adjusted  in  inten¬ 
sity  so  that  all  relaxation  transients  started  with  the  same 
optically  generated  electron  density.  Figure  18  shows  the 
transient  response  of  the  SPM  structures  (similar  results  ap¬ 
ply  to  the  DPM)  as  represented  by  the  electron  density  time 
history  at  a  distance  of  100  A  from  the  center  of  the  precipi¬ 
tate.  The  results  indicate  that  picosecond  lifetimes  are  ob¬ 
tained  for  the  case  where  the  defect  density  in  the  surround¬ 
ing  material  is  at  the  1  X  lO^Vcm^  level.  A  drop  in  the  defect 
level  to  1  X 10 ^^/cm^  increases  the  lifetime  to  tens  of  pico¬ 
seconds.  Further  reduction  in  the  residual  defect  density  has 
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FIG.  19.  Transient  recovery  of  a  Schottky  barrier  precipitate  from  optical 
radiation  as  a  function  of  optical  power  intensity.  Residual  trap  density 
surrounding  the  precipitate  is  10‘*/cm^.  The  solid  line  curve  for  this  figure 
represents  the  same  data  as  the  solid  line  curve  for  Figure  18. 


little  effect  as  the  optically  generated  carriers  diffuse  through 
the  computational  boundary  of  the  simulation  in  approxi¬ 
mately  100  ps.  The  recovery  time  for  the  lO^Vcm^  calcula¬ 
tion  is  very  similar  to  that  of  Figure  7;  indeed  the  depen¬ 
dence  on  trap  density  follows  standard  arguments. 

Additional  enhanced  calculations  were  performed  with 
residual  traps  distributed  as  Saturn-type  rings  around  the  pre¬ 
cipitate,  Enhanced  recovery  of  the  transient  occurred  only  in 
those  regions  where  the  trap  concentration  increased.  Calcu¬ 
lations  were  also  performed  with  different  values  for  the  cap¬ 
ture  coefficients.  In  regions  of  the  shortest  recombination 
times  the  carriers  approached  equilibrium  on  the  shortest 
time  scale.  For  both  the  SPM  and  DPM,  the  time  scales  were 
dominated  by  the  material  properties  surrounding  the  pre¬ 
cipitate,  and  not  the  precipitate.  Some  experiments^^  seem  to 
suggest  that  the  defects  between  the  precipitates  are  negli¬ 
gible.  We  have  no  explanation  for  this  discrepancy. 

For  the  modest  optical  power  used  in  this  study,  the 
lifetime  of  the  transient  photoresponse  is  approximately  in¬ 
dependent  of  optical  power.  Figure  19  shows  this  result  for  a 
residual  defect  density  of  1  X  lO^Vcm^  over  a  range  of  opti¬ 
cal  power. 
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FIG.  20.  Comparison  of  DPM  and  SPM.  Variable  is  the  shallow  donor 
concentration.  Solid  lines  in  the  figure  arc  all  for  the  Schottky  model;  sym¬ 
bols  arc  for  the  defect  model.  Inset  is  for  SPM  with  a  donor  concentration  of 
6xl0**/cm^ 


XVI.  DEPLETION  ZONES,  SPHERICAL 
PRECIPITATES,  SHALLOW  DONORS 

The  results  of  the  study  suggest  that  LT  GaAs  comprised 
of  precipitates  with  defects  in  the  surrounding  material  could 
account  for  much  of  the  dc  and  transient  behavior  of  the 
structure.  Within  this  fiamework  it  is  necessary  to  revisit  the 
role  of  precipitate  spacing  on  the  resistivity  of  the  structure. 
Calculations  were  performed  for  both  the  DPM  and  SPM 
with  a  single  spherical  precipitate  surrounded  by  defects  at  a 
level  of  I  X  lO^Vcm^  with  NyN\=l.  Two  results  were  an¬ 
ticipated:  first,  the  material  remains  semi-insulating  for 
higher  concentrations  of  shallow  donors,  when  compared  to 
the  materials  without  surrounding  traps;  second,  the  deple¬ 
tion  zones  for  the  Schottky  barrier  precipitates  are  broader 
than  those  for  the  defect  precipitates. 

The  results  of  the  simulations  are  shown  in  Figure  20, 
and  should  be  compared  to  the  results  of  Figures  15  and  16, 
which  were  performed  for  spherical  precipitates  without  any 
residual  traps.  The  lines  in  Figure  20  are  the  SPM  results;  the 
symbols  are  for  the  DPM.  As  can  be  seen,  for  distances  far 
from  the  precipitate,  the  residual  defects  control  the  level  of 
free  electrons.  For  the  DPM  structure  at  No  =  5XlO^Vcm^ 


0  500  1000  1500  2000 

Distance  from  barrier  (angstroms) 

FIG.  21.  Conduction  and  valence  band  profiles  under  various  bias  condi¬ 
tions,  using  the  zero-band  gap  semi-metal  model  for  the  Schottky  barrier. 

the  free  electron  concentration  is  only  about  1  X  10*Vcm^;  as 
the  doping  level  is  raised  slightly,  the  electron  concentration 
rises  rapidly.  By  contrast,  for  the  SPM  structure  at  a  donor 
density  of  6X10^’/cm^,  the  free-carrier  density  at  100  A 
from  the  center  of  the  precipitate  is  approximately  2 
X  lO^Vcm^.  Thus,  if  we  were  dealing  with  precipitates  that 
were  separated  by  a  distance  of  200  A  only  the  SPM  would 
yield  semi-insulating  material.  For  this  structure  the  transi¬ 
tion  to  conducting  material  occurs  with  a  donor  density  of 
approximately  6X10*Vcm^.  The  inset  to  Figure  20  shows 
the  density  distribution  for  this  higher  donor  level,  where  at 
100  A  from  the  center  of  the  precipitate  the  density  is  near  its 
maximum  attainable  value. 

XVIL  CONCLUSIONS 

We  have  performed  a  series  of  simulations  aimed  at 
modeling  the  physical  properties  of  annealed  LT  GaAs  and 
gaining  an  understanding  of  the  role  and  character  of  the 
precipitates  which  appear  in  the  material.  The  results  of  the 
study,  as  summarized  in  Table  I,  indicate  that  isolated  calcu¬ 
lations  that  ignore  such  fundamental  details  of  the  material 
characteristics  as  precipitates  can  explain  the  dc  and  pho¬ 
totransients.  Since  the  experimental  observation  of  precipi¬ 
tates  means  they  are  not  to  be  ignored,  when  precipitates  are 
included  in  the  study  it  was  concluded  that  all  of  the  experi¬ 
mental  results  including  transients  can  be  explained,  pro¬ 
vided  the  precipitates  are  surrounded  by  residual  defects.  The 
uniformity  of  the  residual  defect  distribution  is  less  of  an 
issue  than  is  the  concentration  of  residual  defects  and  the 
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FIG.  22.  Electron  density  proSles  under  various  bias  conditions,  corre¬ 
sponding  to  Figure  21. 

values  of  their  associated  capture  coefficients.  While  it  is  not 
possible  to  say  definitively  whether  the  precipitates  in  LT 
GaAs  act  as  buried  Schottky  barriers,  or  whether  the  precipi¬ 
tates  are  clusters  of  defects,  the  experimental  evidence 
strongly  suggests  that  the  precipitates  are  not  electrically 
passive,  and  the  simulations  suggest  that  they  are  more  likely 
to  be  buried  Schottky  bairiers.  The  requirement  of  traps  be¬ 
tween  precipitates  is  not  consistent  with  the  conjecture  of 
Luo  et  as  hopping  conduction  in  annealed  material 

with  intra-precipitate  traps  would  render  the  material  con¬ 
ductive.  Additional  study  in  this  area  is  clearly  called  for. 

By  way  of  demonstration,  the  simulation  results  clearly 
indicate  that  dc  measurements  will  not  adequately  provide  a 
clear  picture  of  the  material  properties  of  LT  GaAs 
— transient  measurements  are  also  required.  The  issue  of 
multidimensional  effects  is  briefly  mentioned  in  the  study.  It 
is  worthwhile  emphasizing  that  after  multidimensional  calcu¬ 
lations  are  performed  some  of  the  conclusions  of  this  article 
may  be  altered,  particularly  with  respect  to  such  features  as 
the  device  length  dependence  of  the  do  IV  measurements. 
For  here  we  expect  that  as  carriers  are  injected  into  the  non- 
stoichiometric  material,  they  will  begin  to  alter  the  size  of 
the  depletion  zones  surrounding  the  precipitate,  and  thereby 
alter  the  voltage  at  which  the  current  reaches  the  trap-filled 
limit.  It  is  not  likely  that  the  arguments  concerning  the  tran¬ 
sients  will  be  altered  by  moving  to  three  dimensions. 
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FIG.  23,  For  the  structure  of  Figure  21,  with  a  nominal  Schottky  barrier 
height  of  0.8  cV,  the  variation  of  the  conduction  and  valence  band  as  the 
defect  density  surrounding  the  barrier  is  changed,  insets  denote  blowups  in 
the  vicinity  of  the  interface  of,  respectively,  the  conduction  band  (for  a 
Fermi  level  of  0.6  eV  below  the  conduction  band)  and  the  valence  band  (for 
a  Fermi  level  greater  that  l.O  eV  below  the  conduction  band). 
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APPENDIX  A:  THE  SCHOTTKY  BARRIER  MODEL 

Prior  to  implementing  the  LT  GaAs  Schottky  barrier  pre¬ 
cipitate  simulation  model,  we  performed  a  number  of  planar 
geometry  simulations  of  a  Schottky  barrier  on  n  and  p  type 
material  to  verify  that  the  model  gave  qualitatively  meaning¬ 
ful  results.  The  buried  Schottky  barriers  were  represented  as: 
(1)  a  zero  band  gap  semiconductor  with  a  specified  offset 
voltage  with  respect  to  GaAs,  in  which  (2)  the  electron  and 
hole  density  are  equal  and  high  and  where  (3)  everywhere 
within  the  material,  except  at  the  region  surrounding  the 
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FIG.  24.  Electron  density  distribution  as  determined  by  the  potential  distri¬ 
bution  of  Figure  23  and  Poisson’s  equation. 


FIG.  25.  Current -voltage  characteristics  for  a  simple  Schottky  barrier  pre¬ 
cipitate,  and  one  surrounded  by  defects  where  £<;(x  =  0)  -  £/. >  1.0  cV.  The 
reference  current  density  is  1.6 X  10^  A/cm^. 


boundary  of  the  zero  band  gap  material  and  the  GaAs  mate¬ 
rial,  there  was  local  charge  neutrality.  We  have  not  tried  to 
fine  tune  the  model  to  duplicate  experimental  results.  (Note: 
the  boundaries  of  the  computational  domain  are  treated  as 
ohmic  contacts,  following  standard  drift  and  difftision  prac¬ 
tices.)  This  appendix  summarizes  the  study. 

The  first  simulations  were  for  a  0.8  eV  bamer  (conduc¬ 
tion  band)  on  n  type  material  doped  to  1  X  lO^Vcm^.  Figure 
21  shows  the  band  structure  for  this  configuration  at  bias 
levels  of  0.0,  0.4,  and  0.9  eV.  The  maximum  depletion  oc¬ 
curs  under  zero  applied  bias,  represented  by  the  dotted  lines. 
The  corresponding  electron  density  is  shown  in  Figure  22. 
These  results  are  similar  to  those  found  in  any  textbook. 

Introducing  a  layer  of  traps  in  the  vicinity  of  the 
Schottky  barrier  alters  the  effective  barrier  height  To  dem¬ 
onstrate  this  a  donor  defect  sheet  density  fixed  at 
1  X  10*^/ cm^  was  placed  across  a  100  A  layer  and  the  accep¬ 
tor  defect  level  was  varied  to  change  the  Fermi  level.  Results 
are  presented  for  zero  bias  in  Figure  23  and  identified  by  the 
position  of  the  Fermi  level.  With  a  trap  ratio  such  that  the 
Fermi  level  was  at  0.8  eV  (Figure  23,  dotted  line),  the  band 
structure  is  effectively  unaltered  from  the  case  of  a  Schottky 
barrier  alone.  When  the  trap  ratio  yields  a  Fermi  level  of  0.6 
eV  (Figure  23,  dashed  line)  or  somewhat  greater  than  1.0  eV 
(Figure  23,  solid  line)  a  local  redistribution  of  charge  occurs 
at  the  interface.  For  a  Fermi  level  of  0.6  eV  there  is  a  narrow 
spike  in  the  conduction  band  (topmost  inset  to  Figure  23, 
where  we  show  conduction  band  energy  over  a  100  A  region 
centered  about  the  interface).  For  a  Fermi  level  >1.0  eV, 
there  is  a  narrow  spike  in  the  valence  band  (bottom  inset  to 


Figure  23).  The  corresponding  electron  density  plots  are 
shown  in  Figure  24.  The  spike  in  electron  density  is  a  Fermi 
level  of  0.6  eV,  and  occurs  within  the  ‘metal’.  For  a  higher 
density  ‘metal’  this  change  in  concentration  would  be  negli¬ 
gible.  For  a  Fermi  level  >1.0  eV,  there  is  a  decrease  in 
density  on  the  ‘metal’  side. 

These  results  indicate  the  sensitivity  of  the  Schottky  bar¬ 
rier  contact  to  the  material  properties  of  the  sunounding  re¬ 
gion.  This  is  demonstrated  further  in  Figure  25  where  we 
display  the  forward  current-voltage  characteristics  for  the 
ordinary  Schottky  barrier  (closed  circles)  and  for  a  barrier 
with  a  trap  layer  which  increases  the  effective  barrier  height 
by  0.2  eV  (closed  squares).  The  solid  lines  are  obtained  from 
the  equation 

log  io(~f/ Reference  Current) 


iogio 


exp- 


[  ^  ^effective  barrier  height] 


nkgT 


where  n  in  the  above  equation  is  the  ideality  facto r.^^  For  the 
calculation  without  traps,  /?=  1,  and  the  extrapolated  effec¬ 
tive  barrier  height  is  0.76  eV,  which  is  somewhat  less  than 
the  input  of  0.8  eV.  For  the  calculation  with  traps  1.09. 
and  the  extrapolated  effective  barrier  height  is  approximately 
0.8  eV.  Notice  that  both  curves  show  current  saturation,  with 
the  trap  layer  structure  leading  into  saturation  at  a  value  of 
bias  approximately  0.2  eV  higher  than  that  of  the  study  with¬ 
out  traps.  The  current  voltage  characteristics  for  the  latter 
also  exhibit  additional  structure  at  the  higher  bias  levels,  that 
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FIG.  26.  As  in  Figure  23,  but  for -type  material.  Inset  is  a  blowup  of  the 
valence  band  for  the  case  when  =0)-£f>  1.0  cV. 


=  M/iO  ■ 


E 

\^0, 


arises  from  the  filling  of  traps  in  the  vicinity  of  the  ‘metal’- 
semiconductor  interface.  Apart  fixjm  this  difference,  both 
curves  show  classical  exponential  behavior  prior  to  satura¬ 
tion. 

Finally,  Figures  26  and  27  show  results  for  the  same 
structures  presented  in  Figures  23  and  24,  but  with  the  semi¬ 
conductor  material  doped  to  I  X  10^^/cm^  p  type.  The  results 
are  as  expected,  and  serve  to  demonstrate  the  applicability  of 
the  model  to  both  n  and  p  type  materials. 


where; 

yu.„o  =  8000  cm^A^s 
/x^o“400  cm^A^s 

i;„s= 8.5X10^  cm/s  (B3) 

i;p5=  l.OX  lO^cm/s 
£o  =  4000  V/cm. 

For  the  case  of  doping  in  n-type  material,  the  constant  mo¬ 
bility  in  the  above  expressions  is  replaced  by: 


APPENDIX  B:  CAPTURE  COEFFICIENTS,  THE 
MOBILIPf  MODEL  AND  PRECIPITATE  DENSITY 

The  baseline  capture  coefficients  used  in  this  study  were: 


8000  cm^/V  s 


(B4) 


C„  =  1 0” ^cm Vs  C  -  1 .3 1  X  1 0 ~ ^cm Vs 

- 2.6 X  1 0"  '^cm Vs  C~  =  1 0 ■  ^emVs.  (B 1 ) 

The  following  field  dependent  mobility  model  is  used  in  this 
study.  For  undoped  material: 


where  the  net  ionized  density  includes  the  shallow  donor  and 
acceptors  and  the  deep  donors  and  acceptors.  Equation  (B4) 
is  a  result  of  fitting  empirically  obtained  data. 

In  calculating  the  precipitate  density  in  the  spherical 
DPM  discussion,  we  assumed  the  following:  precipitate  di¬ 
ameter:  A;  precipitate  density:  /^densif>' 

=  10‘VcmV  defect  density:  /^densm'^  lO'VcmV  Number  of 
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defects/precipitate:  10.  In  a  sphere  of  diam¬ 

eter  60  A,  the  density  of  defects  per  precipitate  is: 
[^dcnsir\' /^density ~8.84X  10*^/cni^. 
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